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Abstract 



j> ■ Cusps of superconducting strings can serve as GRB engines. A powerful 

■ beamed pulse of electromagnetic radiation from a cusp produces a jet of ac- 

m ■ 

. celerated particles, whose propagation is terminated by the shock responsi- 

^> • ble for GRB. A single free parameter, the string scale of symmetry breaking 

7] ~ 10 GeV, together with reasonable assumptions about the magnitude 

6 

plains the GRB rate, duration and fluence, as well as the observed ranges of 



of cosmic magnetic fields and the fraction of volume that they occupy, ex- 



these quantities. The wiggles on the string can drive the short-time structures 

X 

of GRB. This model predicts that GRBs are accompanied by strong bursts 
of gravitational radiation which should be detectable by LIGO, VIRGO and 
LISA detectors. Another prediction is the diffuse X- and gamma-ray radiation 
at 8 MeV - 100 GeV with a spectrum and flux comparable to the observed. 
The weakness of the model is the prediction of too low rate of GRBs from 
galaxies, as compared with observations. This suggests that either the cap- 
ture rate of string loops by galaxies is underestimated in our model, or that 
GRBs from cusps are responsible for only a subset of the observed GRBs not 
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associated with galaxies. 
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I. INTRODUCTION 



Existing models of gamma ray bursts (GRBs) face the problem of explaining the tremen- 
dous energy released by the central engine [[J. In the case of isotropic emission, the total 
energy output should be as high as 1.4 x 10 54 ergs, in case of GRB 990123 with redshift 
z = 1.6. Strongly beamed emission is needed for all known engine models, such as mergers 
and hypernovae, but such extreme beaming is difficult to arrange (see recent discussion by 
Blandford and Rees ||). In this paper we show that emission of pulsed electromagnetic 
radiation from cusps of superconducting cosmic strings naturally solves this problem and 
explains the observational GRB data using only one engine parameter [|J. 

Cosmic strings are linear defects that could be formed at a symmetry breaking phase 
transition in the early universe |J. Strings predicted in most grand unified models respond 
to external electromagnetic fields as thin superconducting wires ||. As they move through 
cosmic magnetic fields, such strings develop electric currents. Oscillating loops of supercon- 
ducting string emit short bursts of highly beamed electromagnetic radiation through small 
string segments, centered at peculiar points on a string, cusps, where velocity reaches speed 
of light. 0,| . 

The idea that GRBs could be produced at cusps of superconducting strings was first 
suggested by Babul, Paczynski and Spergel (BPS) and further explored by Paczynski 
Iflfl . They assumed that the bursts originate at very high redshifts (z ~ 100 — 1000), with 



GRB photons produced either directly or in electromagnetic cascades developing due to 
interaction with the microwave background. This model requires the existence of a strong 
primordial magnetic field to generate the string currents. 

As it stands, the BPS model does not agree with observations. The observed GRB 
redshifts are in the range z < 3, and the observed duration of the bursts (10~ 2 s < r < 
10 3 s) is significantly longer than that predicted by the model. On the theoretical side, our 
understanding of cosmic string evolution and of the GRB generation in relativistic jets have 



considerably evolved since the BPS papers were written. Our goal in this paper is to revive 
the BPS idea taking stock of these recent advances. 

As in the BPS model we shall use the cusp of a superconducting string as the central 
engine in GRB. It provides the tremendous engine energy naturally beamed. Our main 
observation is that putting superconducting cusps in a different enviroment, the magnetized 
plasma at a relatively small redshift z, results in a different mechanism of gamma radiation, 
which leads to a good agreement with GRB observational data. 

GRB radiation in our model arises as follows. Low-frequency electromagnetic radiation 
from a cusp loses its energy by accelerating particles of the plasma to very large Lorentz 
factors. Like the initial electromagnetic pulse, the particles are beamed and give rise to 
a hydrodynamical flow in the surrounding gas, terminated by a shock, as in the standard 
fireball theory of GRB [|TT|] (for a review see 

The string symmetry breaking scale rj will be the only string parameter used in our 
calculations. With reasonable assumptions about the magnitude of cosmic magnetic fields 
and the fraction of volume in the Universe that they occupy, this parameter is sufficient to 
account for all main GRB observational quantities: the duration tgrb, the rate of events 
Nqrb, and the fluence S. 

We begin in the next Section with a brief review of cosmic string properties and evo- 
lution, with an emphasis on the physics of cusps and on the generation and dissipation of 
electric current in superconducting strings. (The discussion of the latter topic in the existing 
literature is often oversimplified and sometimes incorrect, so we review it in more detail than 
we otherwise would.) The GRB characteristics in our model are calculated in Section III, 
and the hydrodynamical aspects of the model are discussed in Section IV. In Section V we 
discuss the diffuse X-ray and 7-ray backgrounds predicted by the model, as well as other 
observational predictions, which include GRB repeaters, bursts of gravitational radiation, 
and ultrahigh-energy particles. 
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II. STRING OVERVIEW 



A. String properties and evolution 

Here we briefly review some aspects of cosmic string properties and evolution, which are 
relevant for the discussion below (for a detailed review and references see ||). 

Strings are characterized by the energy scale of symmetry breaking 77, which is given by 
the expectation value of the corresponding Higgs field, ((f)) = rj. The mass per unit length 
of string is given by 

(i) 

An important dimensionless parameter characterizing the gravitational interactions of strings 
is 

Gfi ~ (77/mp) 2 , (2) 

where G is Newton's constant and mp is the Planck mass. In many models this is the only 
relevant string parameter. 

Numerical simulations of cosmic string evolution indicate that strings evolve in a self- 
similar manner |T2|-|H|]. A horizon-size volume at any time t contains a few long strings 



stretching across the volume and a large number of small closed loops. The typical distance 
between long strings and their characteristic curvature radius are both ~ t, but, in addition, 
the strings have small-scale wiggles of wavelength down to 

/ ~ at, (3) 

with a < 1. The typical length of loops being chopped off the long strings is comparable to 
the scale of the smallest wiggles ©. 

The loops oscillate periodically and lose their energy, mostly by gravitational radiation. 
For a loop of invariant length / the oscillation period is T} = 1/2 and the lifetime is 
77 ~ 1/kgGfi, where k g ~ 50 is a numerical coefficient. 



The exact value of the parameter a in (|3]) is not known. Numerical simulations give 
only an upper bound, a < 10~ 3 , while the analysis of gravitational radiation backreaction 
indicates that a > k g Gji. We shall assume, following HP2] , that a is determined by the 
gravitational backreaction, so that 

a ~ kgGfi. (4) 

Note that in this case the loops decay within about one Hubble time of their formation. 
Then, most of the loops at time t have length I ~ at, and their number density is given by 

m(t) ~ a~V 3 . (5) 

Analysis of string equations of motion reveals that oscillating loops tend to form cusps, 
where for a brief period of time the string reaches the speed very close to the speed of light. 
Near a cusp, the string gets contracted by a large factor, its rest energy being turned into 
kinetic energy. For a string segment of invariant length 51 <C /, the maximum contraction 
factor is ~ 1/51, resulting in a Lorentz factor 

7 ~ 1/51. (6) 

To avoid confusion, we note that cusps were originally defined |16| as points of infinite 
contraction, where the string momentarily reaches the speed of light. Strictly speaking, such 
cusps can be formed only on idealized infinitely thin strings. For realistic strings, the cusp 
development is truncated either by the annihilation of overlapping string segments at the 
tip of the cusp [|TT|-|T9|j or for superconducting strings, by the back reaction of charge carriers 



or of the electromagnetic radiation. However, unless the string current is very large, so that 
the energy of the charge carriers is comparable to that of the string itself, the truncation 
occurs at a very large Lorentz factor and the string exhibits cusp-like behavior. Below we 
shall use the word "cusps" to refer to such ultra-relativistic string segments. 
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Cusps typically form a few times during an oscillation period, but it is possible to con- 
struct (somewhat contrived) loop configurations exhibiting no cusps. Apart from various 
backreaction effects, the motion of loops is strictly periodic, and thus cusps reappear at 
nearly the same locations on the string in each oscillation period. 



Another peculiar feature that one can expect to find on string loops is a kink [20]. It 
is characterized by a sharp bent, where the string direction changes discontinuously. Two 
oppositely moving kinks are produced on a loop at the moment when the loop is disconnected 
from a long string. The kinks then run around the loop at the speed of light. 



B. String superconductivity 

As first shown by Witten ||, strings predicted in a wide class of elementary particle 
models behave as superconducting wires. If some fermions acquire their mass as a result of 
the same symmetry breaking that is responsible for the string formation, then these fermions 
are massive outside the string but are massless inside. If in addition some of these fermions 
are electrically charged, then the strings have massless charge carriers which travel along 
the string at the speed of light. The fermion mass outside the string is m = Xr], where A 
is the Yukawa coupling of the fermion to the Higgs field of the string. Yukawa couplings in 
particle physics models are often very small, so it is not unusual to have m <C 77. String 
superconductivity can also be bosonic, with charge carriers being either spin-zero bosons or 
spin-one gauge particles. Here, we shall consider only fermionic superconductivity. 

An electric field E applied along a superconducting string generates an electric current. 
The Fermi momentum of the charge carriers grows with time as pp = eE, where e is the 
elementary charge, and the number of fermions per unit length, n = pf/2tc, also grows, 
h ~ eE. The resulting current, J ~ en, grows at the rate 

dJ/dt ~ e 2 E. (7) 
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A superconducting loop oscillating in a magnetic field B acts as an ac generator and 
develops an ac current of amplitude 

J ~ e 2 Bl. (8) 

This loop current is not homogeneous; it changes direction along the string and is more 
accurately described as current-charge oscillations. Some portions of the loop develop charge 
densities ~ J . For typical values used in the calculations below, B = 1 ■ 1CT 7 G and 
I = at Q ~ 30 pc, with a = 1 ■ 1CT 8 and to ~ 10 10 yrs the present age of the universe, one 
obtains J ~ 2 ■ 10 5 GeV. 

The local value of the string current can be greatly enhanced in the vicinity of cusps. 
The portion of the string that attains a Lorentz factor 7 is contracted by a factor ~ I/7. 
The charge carrier density, and thus the current, are enhanced by the same factor, so the 
current becomes (in the local rest frame of the string) 

J 7 ~ 7 J - (9) 

The growth of electric current at the cusp is terminated at a critical value J max when the 
energy of charge carriers becomes comparable to that of the string itself, (J/e) 2 ~ \i. This 



gives J max and ~f max as [21 



Jmax ~ 677, ^imax ~ (erj/Jo). (10) 

Alternatively, the cusp development can be terminated by small-scale wiggles on the string 
|22|j . If the wiggles contribute a fraction e C 1 to the total energy of the string, then the 



maximum Lorentz factor is less than (|10f) , and is given by 

Imax ~ e" 1/2 . (11) 

The actual value of 7 max is not important for most of the following discussion. 
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In realistic models, the strings have several fermion species as charge carriers. It can 
be shown that fermions of a given species can move only in a certain direction along the 
string. Thus, the charge carrier species can be divided into left-movers and right-movers. If, 
for example, the applied electric field is directed to the right, it produces positively charged 
right-movers and negatively charged left-movers (and vice versa for the opposite direction 
of E). The left-movers and right-movers usually differ by flavour, lepton number, or some 
other conserved or weakly violated quantum number. 

In the absence of an external electromagnetic field, the current in an oscillating loop 
decays due to various dissipation mechanisms. These include: scattering of left- and right- 
,23|, electromagnetic back-reaction P5|]2E|| , and plasma effects |^7|| . 



movers 



Charge carrier loss due to scattering of left- and right-movers is highly model-dependent. 
If the scattering is mediated by superheavy gauge bosons of mass Mx ~ 10 15 GeV, then the 
characteristic scattering time is p3| 



J 



-5 



r «~ 3 ' 10 llFfevJ < 12) 

For J < 10 4 GeV this time is greater than the age of the universe, but r sc decreases rapidly 
with the growth of the current and becomes comparable to the typical oscillation period 
of loops [T\ ~ 100 yrs for / ~ 30 pc) for J ~ 3 • 10 5 GeV. In near-cusp regions, where 
J 3> 10 5 GeV, charge carrier scattering becomes very efficient. 

We note, however, that this current loss mechanism has an important limitation. The 
densities of left- and right-moving charge carriers are typically not equal, and even if scatter- 
ing were 100% efficient, it would stop after eliminating the minority charge carriers, leaving 
the string with a chiral current (that is, with a current consisting of only left- or right- 
movers). This is what we expect to happen in the vicinity of cusps. 

The electromagnetic back-reaction typically damps the loop current on a timescale r em ~ 
l/e 2 ~ 100/, which is much shorter than the loop's lifetime. It tends to damp the spatial 
component of the current, with the total charge of the loop remaining the same, so in 
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the absence of other effects the end result would be left- and right-moving currents of the 
same magnitude and the same charge.[] Combined with scattering of charge carriers, this 
mechanism can dissipate loop charges and currents, even in the chiral case. Moreover, the 



string charge is almost completely screened by a vacuum condensate fl28| , so the string 
is effectively neutral even if the scattering rate is low and there is some residual charge. 
It should be noted that the physics of electromagnetic back-reaction can be significantly 
modified by plasma effects, which are presently not well understood. Thompson [27] has 
argued that current damping becomes more efficient in the presence of plasma. 

Another mechanism that can dissipate a large chiral current operates when a loop oscil- 
lates in an external magnetic field. The emf induced in the loop oscillates with the same 
period. Suppose for definiteness that the loop initially has a chiral current Jj consisting 
of positively charged right-movers. When the emf is directed oppositely to this current, 
the magnitude of the right-moving current is reduced by ~ Jo an d a positively charged left- 
moving current of magnitude ~ Jo is generated, with Jo from Eq. Left- and right-movers 
can now scatter off the string, and if r sc < TJ, the chiral component of the current will be 
reduced by ~ J . The initial current will then be dissipated in ~ Ji/Jo oscillations. 

The effect of all these dissipation mechanisms is to damp the loop's charge and current 
on a timescale 

T d ~ (1 - 100)/. (13) 

This means in particular that the loop quickly forgets any initial charge or current that it 
inherits when it is chopped off the long string network. The magnitude of the current in a 



1 Spergel et. al. [^6| argued that the dc component of the current cannot be changed by the e-m 
back-reaction. However, their Eq. (11) which they quote in support of this statement is in fact an 
expression of charge conservation. 
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loop is determined mainly by the local magnitude of the cosmic magnetic field, as in Eq. 
®- 

We note finally that Eq. @ for the current is modified when the loop has an appreciable 
center-of-mass velocity v ||. In this case, the amplitude of current-charge oscillations grows 
linearly with time, until the growth is hampered by the damping processes. The resulting 
amplitude is 

J ~ e 2 Bvr d . (14) 

Loops can have high center-of-mass velocities, v ~ 1, but in view of the uncertainty in the 
damping time (|H|) we shall use the estimate (|S|) for the current. 

III. GRB ENGINE 

There are three types of sites in the universe where magneic fields can induce large electric 
currents in the strings. They are compact structures (galaxies and clusters of galaxies), 
voids, and walls (filaments and sheets) of the large-scale structures. The total rate of GRBs 
is dominated by the walls, and further on we shall concentrate on these structures only. 

Magnetic fields in our scenario are assumed to be generated in young galaxies during the 
bright phase of their evolution [^] and then dispersed by galactic winds in the intergalactic 
space. Then at present the fields are concentrated in the filaments and sheets of the large- 
scale structure [[Jl],[n|]. 

Assuming that magnetic fields were generated at some z ~ Zb (galaxy formation epoch) 
and then remained frozen in the extragalactic plasma, we obtain 

B(z) = B (l + z) 2 , (15) 

where the characteristic field strength at the present time B can be estimated from the 
equipartition condition as B ~ 1CT 7 G [j30|1 . 
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With sheets of characteristic size L ~ (20 — 50)/i _1 Mpc and thickness D ~ 5/i _1 Mpc, 
we can estimate the fraction of the space occupied by the walls with magnetized plasma as 
J'b ~ D/L ~ 0.1. For numerical estimates below we shall use zb ~ 4. 

We shall now estimate the physical quantities characterizing GRBs powered by cusps of 
superconducting strings. In what follows we assume that the universe is spatially flat, is 
dominated by non-relativistic matter, and has age t = 0.87 • 10 10 yr, which corresponds to 
dimensionless Hubble constant h = 0.75. 

A. GRB rate and fluence 

Due to the large current, the cusp produces a powerful pulse of electromagnetic radiation. 
The total energy of the pulse is given by 0|| ~ 2k em JoJ max l, where I ~ at is the length 
of the loop, and the coefficient k em ~ 10 is taken from numerical calculations 0. This 
radiation is emitted within a very narrow cone of openening angle 6 min ~ l/^max- The 
angular distribution of radiated energy at larger angles is given by [[7| 

d£ em /dn ~ k em j 2 i/e\ (16) 

For a GRB originating at redshift z and seen at angle 6 with respect to the string velocity 
at the cusp, we have from Eqs.(p|)-(p^|) 

d£ em /dQ ~ k em e 4 aH 3 B 2 (l + z^O' 3 , (17) 

The Lorentz factor of the relevant string segment near the cusp is 7 ~ 1/9. The duration of 
the cusp event as seen by a distant observer is 

t c ~ (1 + z)(at/2) 7 - 3 ~ (a*o/2)(l + z)- 1/2 9 3 . (18) 

One can expect that the observed duration of GRB is tqrb ~ r c- This expectation will be 
justified by the hydrodynamical analysis in Section IV. 
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The fluence, defined as the total energy per unit area of the detector, is |10| 



s~(i + z )(ds em /dn)di 2 {z), 



(19) 



where di(z) = 3to(l + 2 ) 1 ^ 2 [(1 + Z Y^ 2 ~ 1] is the luminosity distance. 

The rate of GRBs originating at cusps in the redshift interval dz and seen at an angle 9 
in the interval d9 is given by 



Here, u(t) ~ ni(t)/Ti ~ 2a 2 t 4 is the number of cusp events per unit spacetime volume, 
Ti ~ at/2 is the oscillation period of a loop, dV = 547r^[(l + z) 1/2 - 1] 2 (1 + z)~ u/2 dz is the 
proper volume between redshifts z and z + dz, and we have used the relation dt = {l + z)dt. 

Since different cusp events originate at different redshifts and are seen at different angles, 
our model automatically gives a distribution of durations and fluences of GRBs. The angle 
9 is related to the Lorentz factor of the relevant portion of the string as 9 ~ I/7, and from 
Eqs.(|TJ),(|lJ) we have 



Here, 70 ~ 190, as = a/10 8 , and the fluence S and the magnetic field B are expressed 
as S = S_ s ■ 10~ 8 erg/cm 2 and B = B_ 7 ■ 10~ 7 G. 

Very large values of 7 ~ jmax, which correspond (for a given redshift) to largest fluences, 
may not be seen at all because the radiation is emitted into a too narrow solid angle and the 
observed rates of these events are too small. The minimum value 7(2; S m i n ) is determined by 
the smallest fluence that is observed, e.g. for GRBs at z ^ 1 with S m i n ~ 3 • 10~ 8 erg/cm 2 , 
7mm ~ 170. Another lower limit on 7, which dominates at small z, follows from the condition 
of compactness U and is given by 7 > 100 (see Section IV). 

The total rate of GRBs with fluence larger than 5* is obtained by integrating Eq . ([20|) 
over 9 from 7mL( 2 to ■~y~ 1 (z; S) and over z from to min[z m ; zb], with z m from 7 max (-2 m ) = 



dN G RB ~ is 




(20) 



7 (z; s) ~ 7o«:^i / 8 3 J B: 7 2/3 [(v / TT^ - i) 2 v / m] 1/3 . 



(21) 
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•y(z m ;S). For relatively small fluences, S_ 8 < S c = 0.03(7 ma:r (0)a:_8/7o) 3 -B 2 7 , z B < z m and 
we obtain 



N GRB (> S) ~ ^ /_ dV(z)(l + zfr\*\ s) 



SB [ ZB ^ rf ^ , ^5„,-2, 

2aH\ Jo 

~ 3 • 10 2 S: 8 2/3 5% 3 yr" 1 . (22) 

Remarkably, this rate in our model does not depend on any string parameters and is deter- 
mined (for a given value of S) almost entirely by the magnetic field Bq. It agrees with the 
observed rate for B_-j ~ 1 (formally, the observed rate N ~ SOOyr^ 1 at S > 1 • 10~ 7 erg/cm 2 
gives -B_ 7 = 3.2). The predicted slope Ngrb(> S) oc S~ 2 / 3 is also in a reasonable agreement 
with the observed one N obs (> S) oc S~ - 55 at relatively small fluences |52| . 



For large fluences SLg > S c , integration of Eq.(]2{J) gives Ngrb{> S) oc S~ 3 / 2 . Observa- 
tionally, the transition to this regime occurs at Ss ~ 10 2 — 10 3 . This can be accounted for 
if the cusp development is terminated by small-scale wiggles with fractional energy in the 
wiggles e ~ 10 _7 a 2 ^B^ . Alternatively, if "f m ax is determined by the back-reaction of the 
charge carriers, Eq. (fL0|) , then the regime (p2|) holds for larger S-s, and observed steepening 
of the distribution at large S can be due to the reduced efficiency of BATSE to detection 
of bursts with large 7. Indeed, large 7 results in a large Lorentz factor 'Jcd of the emitting 
region (see Section IV), and at ■jcd ~ 10 3 photons start to escape from the BATSE range. 

B. GRB duration 

The duration of GRBs originating at redshift z and having fluence S is readily calculated 
from Eqs.flT|) and (gTJ) as 

r GRB re 200 °y 2? (1 + z)-\VT+^ ~ I)' 2 s (23) 



From Eqs. (p0|) and (|18) we find the rate of GRBs with durations in the interval dr and 
redshifts in the interval dz, 
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dN ~ 10 2 a~V ( —) 2/3 (VT+^ ~ 1) 2 (1 + z)- x '«dz—. (24) 

\ Ott J T 

The distribution of GRB durations is found by integrating this over z. The integration is 
restricted by z < z B and S > S min ~ 3 • 10~ 8 erg/cm 2 . The latter condition can be expressed 
as z < z(t), where z(r) is the solution of Eq. ( p3| ) for z with S ~ S m i n . 

The distribution changes its form at the characteristic value r* defined by z(t*) = zb- 
With z B = 4, Eq. © gives 

r, ~ 8.7a 4 _ 8 B 2 _ 7 s. (25) 

For r < t* we have rfiV oc T 2 ^ 3 dr/r, and for r ^> r*, diV oc r~ 5 ^dr/r. We thus see that 
the distribution is peaked at r ~ r*. 

The largest value of r in our model is obtained from Eq. fllSD with 9 ~ 9 max ~ 10~ 2 , 
t grb ~ 10 3 a_g s. There is no sharp lower cutoff, but very small values of r will not be 
observed due to the low rate of events. With the rate ~ 10 2 yr _1 near the peak of the 
distribution, the rate of events with r ~ lO -4 ^ is about 0.1 yr _1 . 

A lower bound on r is also set by the detector resolution (~ 10~ 2 s for BATSE). Hence, 
we have T min ~ max{10~ 4 r^,, 10~ 2 s}. 

The observed distribution of GRB durations extends from ~ 10~ 2 s to ~ 10 3 s. The 
distribution is bimodal, with peaks at 0.5 s and 15 s |33|], and there are some observational 



indications that short and long GRBs may have different origin. Our model is probably 
better suited to describe the short GRB population (see Section V). With t* ~ 0.5 s and 
-B_7 ~ 3, Eq.flSS]) gives as ~ 0.3. This corresponds to the string symmetry breaking scale 
T) ~ 1 • 10 14 GeV. The range of GRB durations is then given by ~ 10~ 2 s, Tq^ b ~ 10 3 s. 

It should be noted that the validity of our simple one-parameter model does not extend 
beyond rough order-of-magnitude estimates (see Section VI). In particular, it is not expected 
to give the correct duration distribution N(t), and identifying the peaks of the theoretical 
and observed distributions may therefore exceed the accuracy of the model. A more con- 
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servative approach is to require that the characteristic duration r* lies within the observed 
range of GRB durations. This gives 0.2 < a_ 8 < 3. 

IV. ACCELERATION AND HYDRODYNAMICS 

A beam of low-frequency e-m radiation propagating in plasma produces a beam of accel- 
erated particles. 

The characteristic frequency of e-m radiation in a pulse produced by a cusp segment with 
Lorentz factor 7 is 

"«» ~ ^7 3 (1 + z? 12 = 4.6( 7 /10 3 ) 3 (l + z f 2 a-_l s~\ (26) 

The plasma frequency in the intergalactic gas of density n = n_5i0 -5 cm -3 , 

Lu p i = 1.8 ■ IOV/5 s _1 , (27) 

is higher than uj em when 7 % 3.4 ■ 10 3 n]/^ a]/g (1 + z)^ 1 ^ 2 . Therefore, the low-frequency 
radiation from the cusp cannot propagate in plasma. In fact, the energy density of e-m 
beam is much larger than that of the plasma, and the beam would push the plasma away 
even in the case u> em > uj v \. This process occurs due to the acceleration of plasma particles. 

Let us consider the propagation of a charged test particle in a strong, low-frequency 
e-m wave. For a time interval much shorter than the period of the wave, t <C l/u em , the 
e-m field of the wave can be approximated by static, orthogonal electric and magnetic fields 



of equal magnitude. Solution of the equations of motion (see e.g. @]) shows that both 
positive and negative charges are accelerated mainly in the direction of wave propagation, 
n = (E x B)/EB, with their Lorentz factor increasing with time as 

^^(w^T- (28) 

where m is the particle's mass. The synchrotron energy loss of an accelerated particle is 
small, because when it moves in the direction of wave propagation, n, the electric force, eE, 
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and the magnetic force ev x B almost exactly compensate each other: e(E + v n x B) « 0. 
This regime of acceleration is practically the same as in the Gunn-Ostriker mechanism [j35f] . 

For an e-m wave in vacuum, a test particle would be accelerated at t ~ \juj em up to 
a very large Lorentz factor. But the maximum Lorentz factor of the beam is saturated at 
the value 75, when the energy of the beam reaches the energy of the original e-m pulse: 
^b m lb ~ 4m> This results in the Lorentz factor of the beam 

lb ~ 4 • 10 2 B 2 _ 7 nZl(l + ^) 4 (7/100) 6 . (29) 

Let us now turn to the hydrodynamical phenomena in which the gamma radiation of 
the burst is actually generated. The beam of accelerated particles pushes the gas with the 
frozen magnetic field ahead of it, producing an external shock in surrounding plasma and 
a reverse shock in the beam material, as in the case of "ordinary" fireball (for a review see 
The difference is that the beam propagates with a very large Lorentz factor 7^ ^> 7, 
where 7 is the Lorentz factor of the cusp (the precise value of 75 is not important for our 
discussion). Another difference is that the beam propagates in a very low-density gas. The 
beam can be regarded as a narrow shell of relativistic particles of width A ~ V^7 3 i n the 
observer's frame. 

The gamma radiation of the burst is produced as synchrotron radiation of electrons 
accelerated by external and reverse shocks. Naively, the duration of synchrotron radiation, 
i.e. tgrb, is determined by the thickness of the shell as tgrb ~ A. This is confirmed by a 
more detailed analysis, as follows. The reverse shock in our case is ultrarelativistic [|36|,[I|1 . 
The neccessary condition for that, pb/p < 7^, is satisfied with a wide margin (here pb is 
the baryon density in the beam and p is the density of unperturbed gas). In this case, the 
shock dynamics and the GRB duration are determined by two hydrodynamical parameters 
[U]. They are the thickness of the shell A and the Sedov length, defined as the distance 
travelled by the shell when the mass of the snow-ploughed gas becomes comparable to the 
initial energy of the beam. The latter is given by ls e d ~ (Siso/p) 1 ^ 3 - 
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The reverse shock enters the shell and, as it propagates there, it strongly decelerates the 
shell. The synchrotron radiation occurs mainly in the shocked regions of the shell and of the 
external plasma. The surface separating these two regions, the contact discontinuity (CD) 
surface, propagates with the same velocity as the shocked plasma, where the GRB radiation 
is produced. 

The synchrotron radiation ceases when the reverse shock reaches the inner boundary of 
the shell. This occurs at a distance ~ Zj^A 1 / 4 when the Lorentz factor of the CD surface 
is 

Icd ~ (WA) 3/8 ~ O.li^Vf (1 + ^/ 2 7 3/2 . (30) 

Note that these values do not depend on the Lorentz factor of the beam 7^ and are determined 
by the cusp Lorentz factor 7. The size of the synchrotron emitting region is of the order 
Ra, and the Lorentz factor of this region is equal to 'Jcd- The compactness condition [|TJ 
requires ■jcd ~ 10 2 , and Eq. ( p0[) yields 7 ^ 10 2 which we used earlier in Section III. The 
duration of GRB is given by 

tgrb ~ #a/27cd ~ V27 3 , (31) 

i.e. it is equal to the duration of the cusp event given by Eq . fll8|) . The energy that goes into 
synchrotron radiation is comparable to the energy of the electromagnetic pulse. 

V. PREDICTIONS AND PROBLEMS 

In this section we shall consider a number of predictions of our model. Some of these 
predictions pose potential problems. 

(i) Short-time structure of GRBs. 

Most of GRBs exhibit a complex short-time structure. These variations must be a prop- 
erty of inner engine [p],p7|. In the cusp model they can be naturally produced by wiggles. 
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Wiggles are amplified in near-cusp regions and, acting like minicusps, produce a sequence of 
successive fireballs. A quantitative analysis of this effect would require a detailed study of 
the gravitational back-reaction, which controls the amplitude of the wiggles. 

(ii) Repeaters. 

Cusps reappear on a loop with a period of loop oscillation, producing nearly identical 
GRBs. In our model, where all loops have the same length I — at at a given cosmological 
epoch t, the recurrence time, TJ ~ (1 + z)at/2 ~ 50a_8(l + z)^ 1 ^ 2 yr, is too long to be 
observed by BATSE and other detectors. In a more realistic models, some fraction of loops 
would have lengths smaller than at and thus shorter recurrence periods. This fraction is 
model-dependent. Moreover, GRBs from repeaters with / < at must be weak and have short 
durations. 

The GRB fluence from a loop of length I produced by a string segment with Lorentz 
factor 7 can be readily calculated as 

After a change of variables from 7 and / to r GRB and the reccurence period T rec = 1(1 + z)/2, 
we obtain 

tore ~ k em e _ [(i + ^ i/2 _ i]2 _. (33) 

A search for repeaters with T rec < 5 yr requires, according Eq. ([3~3"l) , low fluences S ^ 1CT 7 
erg/cm 2 and short durations tgrb ~ 40 ms. The BATSE efficiency is low for such events 
and the repeating burst could easily have been lost. The total number of GRBs shorter 



than 40 ms in BATSE 3B catalogue is less than 5 [3£ 



(in) Host galaxies 

The discovery of GRB aferglows revealed an association of long-duration GRBs with 
galaxies (see 0] for a review). 19 GRBs with long durations are found to be undoubtedly 
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hosted by normal galaxies ||41|| . For 10 of them redshifts are found to be typically 1-3. For 
many bright bursts, which are most probably at small distances, no host galaxies have been 
found. For example, for 16 bright bursts observed by the Interplanetary Network with small 
error boxes, no galaxies are found with magnitudes from 20 to 24. This suggests that some 
of the GRBs are not hosted by galaxies. 

In our model, the fraction of loops captured by galaxes is expected to be small, due to 
the high velocities of the loops. The most straightforward way to reconcile the model with 
observations is to assume that cusps are responsible only for a subset of the observed GRBs 
not associated with galaxies. Such a subset could include the short-duration GRBs, for which 
no host galaxies have yet been detected. With the choice of parameters B^ 7 ~ 3, «_s ~ 0.3, 
as in Section III, the distribution of GRB durations is peaked at tgrb ~ 0.5 s. At the same 
time, the tail of the distribution extends all the way to tqrb ~ 10 3 s, and thus the model 
can account for some of the long GRBs as well. This particular possibility meets another 
problem, since short GRBs do not show deviation from Euclidean distribution. However, 
it is often suggested that GRBs comprise a few subclasses, and the existence of a no-host 
subclass remains plausible. 

An alternative possibility should be also mentined. In string evolution models with 
a ^> k g Gn, the lifetime of the loops is 77 ^> t, so the loops will be slowed down by the 
expansion of the universe and a substantial fraction of them can be captured in galaxies. 

(iv) Bursts of gravitational radiation 

Our model predicts that GRBs should be accompanied by strong bursts of gravitational 
radiation (see also @). The angular distribution of the gravitational wave energy around 
the direction of the cusp is ^| dS g /dQ ~ Gfi 2 l/9, and the dimensionless amplitude of a 
burst of duration r originating at redshift z can be estimated as 

h ~ k-^a^ir/t^il + z)-^z-\ (34) 

or h ~ 10 _21 a% 3 z _1 (r/ls) 1 / 3 for z ^$ 1. Here, we have used the relation F g ~ h 2 /GT 2 ~ 
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(1 + z)(d£ g / dVL) / d 2 L r for the gravitational wave flux and Eq.(fL8D for the burst duration r. 
These gravitational wave bursts are much stronger than expected from more conventional 
sources and should be detectable by the planned LIGO, VIRGO and LISA detectors. It has 
been shown in [42] that gravitational wave bursts from strings are linearly polarized and 
have a characteristic waveform h{t) oc t 1 / 3 . 

(v) X- and 7- ray diffuse radiation 

Tremendous energy [see Eq.([H])] released in a narrow angle 9 ~ l/jmax is not seen in 
GRBs because of the smallness of this angle. The beam of particles accelerated by e-m 
radiation in this narrow cone has a very large Lorentz factor, and the emitted photons 
have energies in excess of 1 TeV. These photons are absorbed in collisions with infrared 
(IR) or microwave photons, collectively denoted as 7^ 7 + 74 — > e + + e~ . Electrons and 
positrons start e-m cascades on microwave photons (755) due to Inverse Compton scattering 
(e + 7bb — > e + 7) and pair production (7 + 7;,;, — > e + + e~). As they degrade in energy, cascade 
electrons are effectively deflected in the extragalactic magnetic field, and the produced diffuse 
gamma radiation is isotropic. The spectrum of remaining cascade photons was calculated 



analytically in [44] (for recent Monte Carlo simulation see [45 1). The analytic spectrum is 
described in terms of three parameters: e 7 , ex and u 1 . 

e 7 is the minimum energy of absorption, i.e. the smallest energy of a photon absorbed on 
IR radiation (e 7 ~ m^/eiR, with the exact value dependent on the spectrum of IR radiation). 
ex is the energy of an IC photon produced by an electron of energy e e = e 7 /2, i.e. by an 
electron born by a photon of energy e 7 . cj 7 is the energy density of cascade radiation. 

The space density of cascade photons, n^(E), is given by flif 



n 7 (E) 



K(E/e x )' 3/2 if E < e x 

K(E/e x y 2 if e x < E < e 7 (35) 
if E > e 7 

where K is a normalization constant which can be expressed in terms of uo cas as 
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The cascade energy density u cas can be calculated as the total energy release in e-m radiation 
of cusps integrated over redshits from up to zb = 4. This gives 

uicas = 5f B k em e 3 r]B /t (37) 

Assuming e 7 ~ 100 GeV due to absorption on IR radiation, we obtain ex ~ 8.1 MeV. Then 
the predicted spectrum in the energy range 8 MeV < E < 100 GeV is 

Itheor(E) ~ 2.5 • lO" 10 ^ £-7 ( 10 3M e v ) cm- 2 s- l sr- l MeV~\ (38) 



This is to be compared with the EGRET flux j|6| for the energy range 5 MeV < E < 
100 GeV, 



E 



-2.1±0.03 



I obs {E) = 1.38 • 10- 9 (^ 10 3 MeV j cm-h-hr^MeV- 1 . (39) 

With a_8 ~ 0.3 and B^-j ~ 3 the predicted flux differs from the observed one by a factor 
of 3. This can be regarded as agreement for an order of magnitude estimate of our simple 
model. 



(vi) Ultra High Energy Cosmic Rays 

GRBs have been suggested as possible sources of the observed ultrahigh-energy cosmic 
rays (UHECR) P?| , ^8| . This idea encounters two difficulties. First, if GRBs are distributed 



uniformly in the universe, UHECR have a classical Greisen-Zatsepin-Kuzmin (GZK) cutoff, 
absent in the observations. Second, the acceleration by an ultrarelativistic shock is possible 
only in the one-loop regime {i.e. due to a single reflection from the shock) jKJ. For a 
standard GRB with a Lorentz factor 7^ ~ 300 it results in the maximum energy E max ~ 
j 2 sh m p ~ 10 14 eV, far too low for UHECR. 

Our model can resolve both of these difficulties, assuming that 7 TOaa; is determined by the 
current backreaction, Eq. (|T0"D . 
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If the magnetic field in the Local Supercluster (LS) is considerably stronger than outside, 
then the cusps in LS are more powerful and the GZK cutoff is less pronounced. 

Cusp segments with large Lorentz factors produce hydro dinamical flows with large 
Lorentz factors, e.g., 7 ~ 2-10 4 corresponds to ■jcd ~ 3-10 5 and E max ~ 'JcD m p ~ 1-10 20 eV. 
Protons with such energies are deflected in the magnetic field of LS and can be observed, 
while protons with much higher energies caused by near-cusp segments with 7 > 10 5 are 
propagating rectilinearly and generally are not seen. A quantitative analysis of the UHECR 
flux in this scenario will be given elsewhere. 



VI. DISCUSSION AND CONCLUSIONS 

The nature of GRB engines is still unknown. There are observational indications that 
they are astrophysical objects: about 25 GRBs are reliably found to be located in galaxies, 
probably in regions of star formation; at least in one case GRB is identified with a supernova 
(SN 1998 bw). The most popular now are astrophysical models, with binary neutron star 
mergers ||50|| , failed supernova [|5l] , hypernova |52 and supranova |k| being the front runners 



(for a critical review see [54]). All these models, however, are not developed enough to give 
quantitative predictions. They also share the difficulty with explaining the large beaming 
factor required for GRBs. 

In contrast, the cosmic string model presented here allows one to obtain quantitative 
predictions for the main observational characteristics of GRBs. In this paper we developed a 
deliberately simplified model, which is characterized by a single free parameter (the energy 
scale 7] of symmetry breaking, or a = kgGrf) and by three other physical quantities, relatively 
well restricted (the magnetic field in filaments and sheets B , the epoch of galaxy formation 
zb, and the density of baryonic matter in the filaments and sheets, a quantity not critical 
for the predictions). Nevertheless, the model correctly accounts for the GRB rate, and for 
the range of GRB fluences and durations. It may also explain the short-time structure of 
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GRBs, the diffuse X- and 7-ray backgrounds, and the ultrahigh-energy cosmic rays. 

The string model predicts recurrence of GRBs with a period of T] ~ 50a_8(l + z)~ 1 / 2 
yrs. Very short bursts may have much shorter recurrence periods, perhaps as short as a 
few years. Observation of these repeaters is a challenge for the future detectors with a high 
efficiency for detection of short bursts. 

Another testable prediction of the model is that GRBs should be accompanied by strong 
bursts of gravitational radiation with a characteristic waveform. 

It must be emphasized that our model involves a number of simplifying assumptions. 
All loops at cosmic time t were assumed to have the same length I ~ at with a ~ k g Gi] 2 , 
while in reality there should be a distribution n(l,t). The evolution law (|T5D for B(z) and 
the assumption of /# = const are also oversimplified. A more realistic model should also 
account for a spatial variation of B. Being basically a one-parameter model, our model may 
predict spurious correlations between the GRB characteristics. In particular, the S oc Tq^ b 
correlation, suggested by Eq.(^3|), holds only at a fixed redshift and tends to be washed out 
when the redshit distribution, the loop length distribution n(l,t), and the inhomogeneous 
spatial distribution of the magnetic field are taken into account. 

Our model meets basically one difficulty: it predicts too low GRB rate from galaxies. 
This discrepancy could be explained if our model strongly underestimates the capture rate of 
string loops by galaxies. For example, if a 3> k g Gfi, then the loops are non-relativistic and 
may be effectively captured by galaxies. Another possibility is that our model could describe 
some subclass of the sources not associated with galaxies. Such a subclass could include the 
short- duration GRBs for which host galaxies are not found. In this case, the model needs a 
smaller a, as discussed in section IIIB. In contrast to the prediction of our model, the short 
bursts do not show strong deviation from the Euclidean distribution. This could be due to 
observational selection effects, since the faint short- duration GRBs which form this subclass 
have a low detection efficiency in BATSE. Alternatively, it could be another subclass of 
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no-host GRBs. 

On the other hand, GRBs from cusps have properties which distinctly distinguish 
them from those produced by collapsars: they are periodically repeating on the scale of 
a few decades for majority of GRBs and on the scale of a few years for faint bursts 
(S ^ 1CT 8 erg/cm 2 ) with short duration tgrb ~ 20 ms. The next generation of GRB 
detectors can examine this prediction. 
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